Diffusion on edges of insulating graphene with intra- and inter-valley scattering 
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Band gap engineering in graphene may open the routes towards transistor devices in which electric 
current can be switched off and on at will. One may, however, ask if a semiconducting band gap 
alone is sufficient to quench the current in graphene. In this paper we demonstrate that despite a 
bulk band gap graphene can still have metallic conductance along the sample edges (provided that 
they are shorter than the localization length). We find this for single-layer graphene with a zigzag- 
type boundary which hosts gapless propagating edge states even in the presence of a bulk band gap. 
By generating inter-valley scattering, sample disorder reduces the edge conductance. However, for 
weak scattering a metallic regime emerges with the diffusive conductance Q = (e 2 /h)(£ KK , /L) per 
spin, where t , is the transport mean-free path due to the inter-valley scattering and L 2> & KK i 
is the edge length. We also take intra-valley scattering by smooth disorder (e.g. by remote ionized 
impurities in the substrate) into account. Albeit contributing to the elastic quasiparticle life-time, 
the intra-valley scattering has no effect on Q. 



T3' 

o. 

<N 

>: 

in 
ov 
t-h ; 

00' 

o: 

(N 



X 



PACS numbers: 



I. INTRODUCTION 



Unique electronic properties of single-atomic-layer 
graphene^ stem from its two-dimensional (2D) semi- 
metallic energy spectrum with gapless conical conduction 
and valence bands. Engineering of a semiconductor- type 
band gap is expected to provide another desirable means 
of controlling electric current in graphene, laying the ba- 
sis for electronic applications 2 -. Several mechanisms of 
the gap generation in the single-layer graphene have been 
discussed in literature 2 - - — , including breaking the sublat- 
tice symmetry on a hexagonal boron-nitride substrate 
(see e.g. Ref. 0) and using mechanical strain (see e.g. 
Refs. 

Although the opening of the band gap could have a de- 
sirable effect on transport in the 2D bulk of the material, 
it is, generally, not sufficient to control electric conduc- 
tion near sample boundaries. Boundaries of graphene 
are natural extended defects that can host unusual elec- 
tronic states such as edge states appearing on a zigzag 
boundary 8 . Various manifestations of such edge states in 
transport properties and spectroscopy of graphene have 
been discussed in recent years (e.g. Refs. l9422|2J-|37j) . 
Most essential for our present discussion is the finding 
that the edge states remain conducting even when the 
2D bulk turns into a band insulator, e.g. due to a stag- 
gered potential breaking the sublattice symmetry (see 
e.g. Refs. l3fj||34l ). This is indeed expected since the 
edge states reside on one of the graphene sublattices, 
and, therefore, the influence of the staggered potential is 
reduced to an energy shift without dramatic changes in 
the edge-state dispersion. Thus, the edge states present 
a potential obstacle for the realization of the band insu- 
lator regime in graphene, providing pathways for leakage 
current. It is of both theoretical and practical interest 
to identify the factors that may help to reduce the edge 



conduction. As one of such factors, in this paper we 
theoretically consider structural disorder involving both 
smooth potential fluctuations, which couple states within 
the same graphene valley (intra-valley scattering), and 
atomically sharp defects generating inter-valley scatter- 
ing. 

Earlier, the influence of disorder on the edge trans- 
port was studied numerically for the conventional semi- 
metallic state of graphene (see e.g. Refs. [9 U2fj| - |23h . Inter- 
estingly, the edge transport cannot be quenched by usual 
potential disorder, e.g. by smooth potential fluctuations 
due to remote ionized impurities. Only atomically sharp 
defects suppress the edge transport by mixing counter- 
propagating edge channels through inter-valley scatter- 
ing. In this paper we extend these findings to the edge 
transport in the insulating graphene with the staggered 
potential, A. Instead of using numerical approaches, we 
perform analytic diagrammatic calculations of the edge 
conductance, explicitly proving that the edge-transport 
mean-free path £kk' is limited only by scattering be- 
tween graphene's two valleys, K and K' ', and unaffected 
by smooth potential disorder. Our calculations indicate 
that the undoped insulating graphene can host diffusive 
metallic edge states with the conductance given per spin 
by 



g = 



h L 
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n v 



Aw(k^ 



k±=± — , (1) 
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where the transport mean-free path £kk' is expressed in 
terms of the Fourier transform of the intervalley-disorder 
correlation function, w(k), and the edge-state velocity v, 
and k± are the edge-state Fermi points relative to the K 
and K' valleys, respectively [L is the distance between 
the source and drain, assumed much larger than £kk', 
but smaller than the typical size of the system that ex- 
hibits full localization 23 ]. The specifics of the insulating 
graphene lies in the fact that the Fermi points k± are 
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shifted with respect to K and K' points by the stag- 
gered potential A. Given furthermore the low edge-state 
velocity metallic diffusive regime under weak 

scattering condition \k±\£ KK , ^> 1 emerges without any 
doping of the material. This is in stark contrast with the 
conventional metallic transport which occurs when the 
Fermi level is pushed into a conduction or valence band. 

The subsequent sections give a complete account of our 
theoretical approach: In Sec. Owe introduce the model 
for the edge states in disorder-free insulating graphene 
with the staggered potential and calculate the edge-state 
Green's functions. In Sec. IIIII we introduce the model of 
disorder, calculate the disorder-averaged Green's func- 
tions, the renormalized edge velocity and, finally, the 
edge conductance from Kubo formula. Section IIVI sum- 
marizes our results. 



II. EDGE STATES IN INSULATING 
SINGLE-LAYER GRAPHENE 

A. Boundary problem 

We begin by analyzing the edge states in disorder-free 
graphene described by the effective four-band Hamilto- 
nian: 



H = v t z erp + At z a z 



(-ihd x ,-iHd y ,0), (2) 



where Pauli matrices r z and u x ,y,z represent the valley 
and sublattice degrees of freedom, respectively [through- 
out the paper products of r - and a- matrices should 
be understood as direct products], vq is the bulk Fermi 
velocity determined by the nearest-neighbor-hopping en- 
ergy and lattice constant, and A is the staggered (e.g. 
substrate-induced^) sublattice potential. Equation @ 
adopts the following convention for the basis states: 



(3) 



where A, B and ± label the sublattices and valleys, re- 
spectively. 

Following our previous work on the edge states in semi- 
metallic graphen o 15 ' 25 ' 27 ' 29 we will work with the Green's 
function, g(r, r'), defined by the equation 




[ei-H]g(r,r')=id(r-r'), 



(4) 



where energy e includes an infinitesimal imaginary part 
iS, with S > (S < 0) for the retarded (advanced) 
Green's function, and I — tqctq is a unit matrix com- 
posed of the unit matrices in valley (to) and sublattice 
(<7o) spaces. Equation (@| will be solved in a semispace 
— oo < x < oo, < y < oo with a single edge at y = 



described by the boundary conditio n 16 i 38 (see also Ap- 
pendix [A}: 



TO 



9\y=o 



TO 



-ern 



i. 



s|j/=o- (5) 
(6) 



It involves two unit vectors n± , orthogonal to each other 
and to the vector normal to the boundary ne \ \y, en- 
suring the vanishing of the particle current normal to 
the edg o 16 ' 38 . The vector components n x or n z serve to 
parametrize the boundary types considered below in Sec. 

urci and irrpi 

The only restrictions on the boundary condition (J5]) 
are the time-reversal symmetry and the absence of the 
intervalley coupling. Thus, Eq. (|S|) can be seen as a 
generalized continuum model for zigzag-type edges. It 
cannot be applied to the armchair edges because the lat- 
ter couples the K and K 1 valleys. However, extended 
armchair edges are unlikely to occur because they are 
less stable than zigzag edges (see, e.g., Ref. H^). As to 
the short-length armchair edges, they can be treated as a 
special type of the boundary defects causing inter-valley 
scattering, which is considered later in Sec. IIIII Therefore, 
Eq. is a good starting point for analyzing defect-free 
graphene. It is also easy to verify that Eq. ((5]) ensures 
vanishing of the normal component of the particle current 
j y (x,0). 



B. Green's function solution 

The Green's function is block-diagonal in valley space, 

(7) 



.9 



g+ 

o g- 



where ^ are 2x2 matrices which, in the basis defined 
by Eq. ([3]), have the following structures: 



9ba 9bb 



g- = ( g BB 9ba ) _ (8 ) 
9ab 9aa 



Let us first calculate the matrix elements of g + . Expand- 
ing 

g+(r,r') = '£g+(y,y>)e ik ^-^/L, 



where L is the edge length, and writing Eq. (j4]) in com- 
ponents, it is straightforward to express the off-diagonal 
elements g\ B \ k and g^ A \ k in terms of the diagonal ones 
as follows 

+ D P- „+ 



9k — v qp+ „+' 

\ e+A 9AA\k 



9AA\k I=K»BB\k \ p± = Hk±hd y , (9) 
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where g^A\k an< ^ 9~BB\ k satisfy the equations: 

e + A 



l d l-i 2 }9AA\ k (y>y') - -j^j 



K - q 2 ]gt B M) 



o 

e-A 



Ik 2 



A 2 - e 2 
h 2 v 2 



s(y-y'), (io) 
S(y-y'), (ii) 



(12) 



The boundary conditions for Eqs. (fTCfl) and (fTTj) follow 
from Eqs. ©, (J7J and © as 



9 

^jtA^ ^') = KA0AA| fc (O,I/'), 

1 - n z e + A 
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(13) 
(14) 



j^9BB\k^y') = K B9B Blk (0,y'), (15) 



l + n 2 e-A 
k b = h k. 

n x hv 



(16) 



We seek the solutions to Eqs. (fTO)) and (ITT1) in the form 
of the linear combinations: 



9BB\ k (y>y') 



A 



2fi 2 v 2 9 
e-A 
"2F^ 



e -«|v-s + CA(y')e~ qv 
e -i\v-v'\ +C B (y , )e- qy 



where the first terms are the Green's functions of the un- 
terminated sublattices, while the second ones are the de- 
caying solutions of the corresponding homogeneous equa- 
tions. The coefficients Ca,b are obtained from boundary 
conditions (fl~3|) and (fTST) with the following results: 



A 



e -q(y+y ) _ e -«|»-» I 



9aA\M)~ 2h 2 v 2 q 

(l + n z )(q + k) + n x (A + e)/hv q ( y+y ') 



2(e — n z A + n x hvok) 
e-A 



(17) 



i{v+y 



-q\y-y I 



0flB|fe(V>tf') - 2h 2 V 2 q 

(1 - n^Xg - fc) + n x (A - e)/fao e _ 9 (j,+„') 
2(e — rt 2 A + n x hvok) 



(18) 



The first terms in Eqs. (|17[) and (| 18[) vanish at the bound- 
ary y — and do not have poles within the gap, |e| < A, 
implying that the edge states are entirely described by 
the second terms. The latter have a pole within the gap 
at e = n z A — n x hv k. Assuming that the energy is close 
to this pole, we can neglect the first terms in Eqs. (fTT)) 
and (|18p and find a compact expression for the matrix 
g+ [see Eq. ©]: 

- + / t\ , Q(n z k + n x A/hv Q ) -aC v+y ') 

9k (y>V) = (0O-0-I1+ — — 7—rqe wry). 

e — n z A + n x nvok 

(19) 




(b) Mass confinement 




FIG. 1: Edge states £+(k) and £-(fe) within bulk band gap 
[—A, A] in valleys K and K' , respectively: (a) zigzag bound- 
ary [see Eqs. ([24]) and ([25]) with n z = 0.99] and (b) mass 
confinement [see e.g. Eq. (|29[) ]. Edge momentum k is mea- 
sured in units of A/hvo. 



The expression for g k can be obtained from Eq. (|19[) by 
replacing vq — > —vq, A —A, and n z — > —n z : 

i\ i <d{-n z k + n x A/hv ) - g ( v + v >) 

g k {y,y) = {ov-(m-) T—rl'i qu+y> . 

e — n z A — n x hvQk 

(20) 

The poles in Eqs. (TT9"]) and (|2T)]) yield the edge spectrum in 
valleys K and K' under the condition that the arguments 
of the Heaviside (theta) functions in the numerators are 
positive. These restrictions are, in turn, enforced by the 
positiveness of the inverse decay length, q. Consequently, 
the edge spectrum is given by 



;+ fc = n z A — n x hv k, n z k + n x - — > 0, for valley K, 

hv 



(21) 



A 



9k = I 9k(y>y) d y 



£_ fc = n z A + n x hvok, —n z k + n x - — > 0, for valley K' . 

hv 

(22) 

These edge states can also be viewed as the solution of a 
strictly ID problem described by Green's functions f| 19[) 
and (|20[) integrated over the transverse coordinate y: 

<Jq - ern± Q(±n z k + n x A/hv ) 
2 e — n z A ± n x hvok 

Below we analyze these results for two commonly con- 
sidered confinement types, viz. the zigzag edge and the 
mass confinement. 



C. Zigzag edge 

In our parametrization of the boundary condition, the 
zigzag-type edge corresponds to the limit n z — > 1. In 
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this case the edge-state spectrum, Eqs. 
reduces to 



e±(k) = A=F hvk, ±k>0, 
where v is the edge-state velocity, 



IB and (122 



= VoV 1 - n l < v o- 



(24) 



(25) 



It vanishes for n z = 1, which corresponds to the flat 
edge-state band. In what follows we keep the velocity v 
finite (but small compared with the bulk velocity vq), so 
that the edge states remain dispersive. We also note that 
the staggered sublattice potential shifts the edge states 
by energy A such that they cross the mid-gap energy at 
finite wave- vectors k± [defined in Eq. (JTJ] with respect to 
the K and K' points (Fig. QJi). The energy shift reflects 
the fact that the zigzag edge states reside on one of the 
sublattices. This is seen from the matrix structure of 
the edge Green's function [Eqs. ([7} and (j2"3")) ] which for 
n z — > 1 has only two diagonal nonzero matrix elements. 
To simplify the model, from now on we will work with 
the effective ID Green's function^ 



9k 





V g- k J 



P+P+g+ + p-p-g^, (26) 



where P^ 1 and P^- are the projector matrices in the valley 
and sublattice spaces, respectively: 



TO ± T z 



PI 



matrix elements g k and g k are given by 

e(±fc) 



e - A ± hvk ' 



(27) 



(28) 



Equation (|2"5j) describes the left-moving (+) and right- 
moving (— ) states originating from valleys K and K' , 
respectively. Comparing Eqs. (J7J and flSJ) with (J25]), we 
see that both left- and right-movers reside on the same 
sublattice A (the case of sublattice B corresponds to the 
boundary condition with n z — > —1), as we mentioned in 
the introduction. 



D. Mass confinement 

This confinement type^ is realized in the limit n z — > 0, 
resulting in the edge-state spectrum, 



(29) 



These edge states have the velocity equal to the bulk 
one, vq, and cross the mid-gap energy at points K and 
K' . Unlike the zigzag edge, there is no energy shift due 




FIG. 2: Diagrammatic representations of (a) Kubo for- 
mula (|32|l for the conductance, (b) Dyson equation (I36[) for 
disorder-averaged Green's functions (thick lines) and (c) equa- 
tion (|45|l for renormalized current vertex in ladder approxima- 
tion. Thin and dashed lines correspond to the unperturbed 
Green's function and disorder correlator, respectively. 



to the staggered potential because in this case the edge 
states propagate on two sublattices. This is again seen 
from the matrix structure of the edge Green's function 
[Eqs. and ([23])] which for n z — > has both diagonal 
and off-diagonal matrix elements in each valley: 



9k 



( 



V 



9k 

-at 
o 





-9t 
9t 










9k_ 
-9 k 






-9k 
9 k 



J 



P+(cr - &x)g~k + P T (CO - <7x)9. 



k • 



(30) 

EH) 



where the valley projectors i~ are defined in Eq. 
and matrix elements g k are given by 

Equations ((30]) and |3l|) recover the corresponding results 
of Ref. 40 for the edge states in 2D topological insulators. 
Their transport properties are reviewed in detail in Ref. 

SI 



III. EDGE CONDUCTANCE 

In the rest of the paper we focus specifically on trans- 
port properties of the edge states on a zigzag-type bound- 
ary. In the metallic regime the conductance of the edge 
states is given by Kubo formula: 



g = 



e 2 h 
2nL 



dk 
2tt 



vG«ie = 0)V k G k A (e = 0) , (32) 
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where G k ' A are the disorder-averaged retarded and ad- 
vanced edge Green's functions, v and V k are the bare and 
renormalized edge-state velocity operators, and L is the 
length of the edge (e.g. the distance between the source 
and drain contacts). To evaluate Eq. (1521) we need to 
calculate first G k ' A and V k in the presence of disorder. 



A. Potential and inter-valley coupling disorder 

We assume that edge disorder can be described by 
equation 

V(x) = IV(x)+ E ( 33 ) 

fi—x,y a— 0,x,y 

where the first diagonal term accounts for smooth ran- 
dom potential fluctuations, e.g. due to remote ionized 
impurities in the substrate. The potential V(x) is char- 
acterized by the correlation function 



(V(x)V(x')) =u(x-x'), 



(34) 



with averaging (...) over disorder configurations (e.g. over 
impurity coordinates). The second term in Eq. (|33[) de- 
scribes atomically sharp defects (lattice vacancies, short- 
length armchair edges etc.) which couple the valleys 
(and, generally, the sublattices) . For this disorder type 
we use the correlation function 



{W lia {x)W lx ,' a '{x')) = w(x - x^S^S, 



fifi 1 u aa' i 



(35) 



where the Kronecker symbols imply completely uncorre- 
cted valley r M and sublattice a a disorder components. 
The calculations below are done for spatially isotropic 
disorder with u(x — x') = u{x' — x) and w(x — x') = 
w(x' — x). 



with two unknown scalar functions G k and G k . The 
ansatz (J3SJ) is valid only for the zigzag-type edge [Cf. 
Eq. ([26]) ] . Inserting this into Eq. (|37| we have 



St. — 



— [ Uk _ k ,Gt,+ ]T w k _ k ,G- k ,] 



r:r; I ^{ Uk _ k ,G k , + ]T w k _ k ,G+] 



dk' 



dk' 



(39) 



Here the first two terms include the potential scatter- 
ing and inter-valley coupling between the left and right 
movers, G k , and G k , . This coupling originates from those 
disorder terms which swap both the valleys and sublat- 
tices. The other two terms in Eq. (131))) result from the 
disorder which swaps either the valleys or the sublattices. 
The latter has no effect on G k because upon inserting 
Eq. (131?)) into Eq. (|5rJ|) the corresponding products of pro- 
jectors vanish: P^P~^P~P+ = 0. Notice that the sum- 
mation over the valley and sublattice indices in Eq. (|3"9")l 
yields the factor of 4. 

Inserting Eq. (|3U)) into Eq. (|36p and collecting coeffi- 
cients at P^Pf we obtain algebraic equations for G k : 



G + k = 



1 



G7 = 



1 



hv(k — k + ) — ' k Hv(k- — k) — S fe 



,(40) 



where are scalar functions given by 



dk 1 



(41) 



B. Disorder-averaged Green's function 



Within the standard self-consistent Born approxima- 
tion^ the disorder-averaged Green's functions can be ob- 
tained from the Dyson equation (see diagram in Fig. ^p) 
which (with suppressed superstripts R, A for brevity) is 
given by 



G k — g k + g k T, k G k 



(36) 



where T, k is the self-energy: 



% = / ^— [uk-k'Gy + Wk-k' ^T^aaGk'T^aa], (37) 

and Uk-k' an d Wk—k' are the Fourier transforms of the 
correlation functions u[x — x') and w(x — x') [see Eqs. 
pi)) and pg])]. We seek the solution to Eq. in the 
form of the projector expansion: 



G k =P?P?Gt+p-p-G k , 



(38) 



Like in conventional metals (see e.g. Ref. |42| ). Eqs. (|4"0")) 
can now be solved using the sharpness of the Green's 
functions G k near Fermi points k±, which yields the fi- 
nite quasiparticle life-time r: 



1 



hv(k - fc + ) + §£ 



IT>( 42 ) 



where the sign s is positive or negative for the retarded 
or advanced functions, respectively. The potential and 
inter-valley scattering mechanisms give additive contri- 
butions to the spectral broadening: 



1 



1 



1 



KK' 

2-kN 



Uq- 



2nN 

4w k k _ — — , 

h 



(43) 
(44) 



where N is the density of states per valley and spin. Since 
potential disorder cannot cause backscattering, the cor- 
responding scattering rate 1/t v involves the correlation 
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function Uk-w at zero momentum transfer k — k' = 0, 
which corresponds to forwardscattering. In contrast, 
inter-valley scattering occurs between the Fermi points 
k± (see Fig. [1]), so that the corresponding scattering rate 
1/t , involves the correlation function Wk + -k- with fi- 
nite momentum transfer k+ — fc_ = 2A/hv. 



The required sum of the renormalized velocities at Fermi 
points k = k± is obtained from Eq. (f50| as 



(52) 



VZ + V,, = v = v 



yielding, finally, the edge conductance: 



C. Vertex renormalization 

We demonstrate below the interplay of the scattering 
times t, t v and t kk , [Eqs. P5|) and (pH)) ] in the disorder- 

renormalized velocity Vk, which is one of the central re- 
sult of this paper. In order to calculate the renormalized 
velocity Vk in the conductance formula (1521 we consider 
the vertex equation in the usual ladder approximation 
(see e.g. Ref. |42j and diagram in Fig. [2};): 

/fik' 
— [uk-k'G^G^ (45) 

+ Wk-k' y^^T^a a G^,Vk'G^,T^(T a }, 



fi.a 



V 



v (p+p+ _ p-p- 



(46) 



Like the edge Green's function (|2"B1) the bare edge velocity 
matrix v (I46[) has only two diagonal elements v and — v 
corresponding to two counter-prapagating channels from 
different valleys. We seek the solution to Eq. (l45l in the 
form of the projector expansion: 



i (47) 

with four unknown scalar functions Vt, VZ, V^ and 

V| 2) . Inserting Eqs. pg j) . (|35 )l and ([37 ]) into Eq. (05) and 
collecting the coefficients at the projectors, we find that 
V^ 1 ' 2 ' can be expressed through V^ by means of 



>+ p-i;( 2 ) 



V, 



(1,2) 



dk' 

~2^ 



2 Wk ^G±«V±G± A , 



(48) 



and V^ satisfy the following equations: 



(49) 

Here again the integral over k' can be calculated using 
the sharpness of the Green's functions at k = k± [see 
Eqs. (|42|) - (|43|) ]. which yields the following result: 



e 2 T KK ,vT e 2 vr^ 



IV. RESULTS AND CONCLUSIONS 

We have demonstrated that in single-layer graphene 
with insulating bulk zigzag-like edges provide pathways 
for metallic conduction. Both intra- and inter- valley scat- 
tering have been taken into account. Although both 
scattering mechanisms contribute to the elastic broad- 
ening of the spectrum P5|) . only the intervalley scat- 
tering time t kk , (|44|) enters the edge conductance ([53]) . 
The transport mean-free path can then be identified as 
I , = vt . We emphasize that Eq. (I5T1) holds under 
weak scattering condition \k±\£ KK , = AT KK ,/h ^> 1, en- 
forced by the staggered potential A in zigzag-type termi- 
nated graphene. Let us discuss qualitatively the depen- 
dence of the edge conductance on the staggered potential 
A and disorder strength, assuming a Gaussian correlation 
function for the inter- valley disorder, 



w(k) = V2irW 2 L c exp(-k 2 L 2 c /2), 



(54) 



where L c is the disorder correlation length and W is the 
root-mean-square amplitude of the disorder. From Eqs. 
CP) and (HI) we have 



0(A) = 



exp 

h 4V2^LL C W 2 V h v 



(2A 2 L 2 



(55) 



We see that the edge conductance exponentially increases 
with A. The reason is that the inter- valley scattering in- 
volves the finite momentum transfer k + — fc_ = 2A/hv 
between the Fermi points k± , with the scattering proba- 
bility cx w(k + — fc_) = w(2A/hv) reducing with A. On 
the other hand, Q gets suppressed algebraically as 1/W 2 
with increasing root-mean-square disorder amplitude W, 
which may help in practice to reduce the edge conduc- 
tance. 



V fe ± = « 



Uk-k± 



u Q + 4w k+ -k- 



-V, 



Awk- 



k = 



u Q + 4w fc+ _ fc „ 



In parallel, we perform the fc-integration in the conduc- 
tance formula ([32]) . obtaining 
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Appendix A: Derivation of boundary condition ([5j) 

We begin by deriving the boundary condition ([5]) for 
the Green's function at the sample edge y — 0. To be 
concrete we consider the retarded Green's function in real 
space and time g(rt,r't'). It is a 4 x 4 matrix in valley and 
sublattice space, with matrix elements g^y (rt, r't') stan- 
dardly expressed through the annihilation tyj(rt) and 
creation ^t,(r'i') field operators as 



x«* i (rt)4(rY) + 4(rY)%(rt))) ! (Al) 

where j (and independently j') runs over the index set 
A + , A_ , £?_ and B + of the basis states introduced in Eq. 
([3]) , the double brackets denote averaging with the equi- 
librium statistical operator, and @(t) is the Heaviside 
function. 

It is obvious from Eq. (|A1[) that the boundary con- 
dition for the Green's function is just the same as for 
^j(rt). Indeed, the boundary condition for *&j{rt) (de- 
rived earlier in Refs. fl6H38l ) can be written as 

* J -(rt)| y =o = M ii * i (rt)| B=0 , (A2) 

where Mji are the elements of the 4x4 matrix 

(If T ° + Tz T -T z 

M = <rn + ern , (A3) 

which is defined in Eq. ([5]) and main text. Inserting Eq. 
(|A2|) into Eq. (|A1|) at y = 0, we obtain the boundary 



condition for the matrix elements of the Green's function: 
g jf (rt,r't') y = = M Jl g ir {rt,r , t% =Q , (A4) 
which along with Eq. (IA3|) yields Eq. (JSJ) - 



Appendix B: Derivation of boundary conditions 
and (TT5|) 

In order to derive these equations we start with the 
boundary condition for the upper block of the Green's 
function, g + \ y= a = — crn + g + \y = o [see Eqs. , and 
©], which has explicit form 

( 9aa 9ab\ = ( n z -n* \ ( 9aa 9 + ab \ (B1 
V 9ba 9bb ) y =o \ ~ n * - n * ) \ 9b a 9bb J y J 

Therefore, for the upper diagonal element we have 

9aa\v=o = n zgA A \y=o - n x9BA\y=o- ( B2 ) 

Expanding in plane waves e ik ( x ~ x ) and using the relation 
[seeEq. ©], 

9$ A \ h M) = ^ K (k + dy)g^ Alk (y,y% (B3) 
we obtain from Eqs. (|B2[) and (|B3|) a closed boundary 
condition for 9jiA\k\v— Q i wn i cn after elementary algebra 
yields the boundary condition (pT3|) . Repeating step by 
step the same calculation for gg B \ y =o in Eq. (|B1|) leads 
to the boundary condition (fl5|) for the other diagonal 
matrix element of the Green's function. 
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